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I. INTRODUCTION 


The origin of tiny but no-zero neutrino masses observed by neutrino oscillation experiments |l|] re¬ 
mains so far a mystery, and so provides us an opportunity to search the new physics beyond the stan¬ 
dard model(SM). Perhaps the simplest scenario which may understand the neutrino puzzle is to introduce 
the Majorana mass that breaks the global B — L symmetry though the dimension-5 Weinberg operator 
ALL$$/a| 1 This effective operator can be realized though various pathways which depends on the 
new physics scales the associated massages lie at. For instance in the case of widely known type-I seesaw 
mechanism A] with right-handed Majorana ingredients Nr as massages, one needs the super-heavy masses 
for Nr i.e 10 1 l ~ 16 GeV to fit the observed sub-eV neutrino mass. The right-handed neutrinos are too heavy 
to be detected at future experiments. In contrast, in the so-called low-scale scenarios, the small neutrino 
mass is not only due to the massage state with heavy mass but also to another naturally small mass param¬ 
eter which breaks the lepton number symmetry This is the basic ideas behind many schemes including 
type-II seesaw model]dj], inverse seesaw model[5 ] and linear seesaw model|6, 7]. In these models, the mass 
of massage particles can be lowered down to TeV or even hundreds GeV, a scale to be explored at collider 
experiments. 

On the other hand, the Planck data has shown that 26% of the energy density of our universe is occupied 
by dark matter. In the view of particle physics, the weakly interacting massive particles(WIMPs) are the 
most promising dark matter candidates. In recent years, a class of models are proposed to incorporate 
the neutrino mass puzzle and the existence of dark matter in a unified framework. In these models, the 
neutrino masses are generated at loop level and the dark matter is naturally contained as a inert particle, 
where the Z 2 symmetry or a U(l)x symmetry is used to guarantee the stability of dark matter. The radiated 


generation of neutrino mass has been realized at one-loop level[8-10], two-loop level[llj-15] and three- 


loop level[16]. T 


-10 


le systematic analysis of one and two-loop realization for with possible topologies are 


performed in Ref] 17]. 


In this paper, we proposed a radiated linear seesaw model where the lepton number violation is due to 
the spontaneous symmetry breaking(SSB) of U(1)r-l gauge symmetry, while the naturally small mass 
parameter is generated at one-loop level. The linear seesaw model was fist studied in the left-right theory 
with gauge group SU (3) x SU (2% x SU (2 )r x 17(1 )b-z,]^], and subsequently inspired on 50(10) theory 
in the presence of gauge singlets (TO. In linear seesaw scenario, T r and T'z are added onto SM so that the 
Lagrangian is given by 


L = MrVl^Ir + M^r^l + + h.c 


( 1 ) 
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The neutrino mass matrix in the basis of (vl, x l ,r lr T l) is 


( 0 

m d 



Ml 

0 

Mqi 

(2) 

\ t*L 

Ml 

0 J 



The light neutrino mass is given by rn v ~ fi i Mp 1 . Note that the ///_ violates the lepton number sym¬ 

metry and plays the role of a naturally small parameter. Thus it seems natural that there exists a suppression 
mechanism where the /i / term is generated via loop diagram meanwhile the soft-breaking of lepton number 
symmetry may attribute to the SSB of B — L gauge symmetry. Moreover, when the WIMPs as the dark 
matter candidates are involved in the loop diagram, we can reasonably assume they are generated by a vac¬ 
uum expectation value with respect to the SSB of B — L gauge symmetry at TeV scale. These are the main 
motivations of this work. 


Following the spirit of Ref. [ 14[], the value of B — L charges should be carefully assigned since the 


anomalies cancelation must be satisfied. It is founded that some new particles may have exotic values of 
B — L charge such that there exists residual Z 2 x Z' 2 symmetry even after SSB of B — L gauge symmetry. 
The Z 2 x Z!) discrete symmetry stabilizes the these particles from decaying to SM ingredients. Thus the 
lightest particle with the same exotic value of B — L charge can be a dark matter candidate. In practise, 
we introduce two classes of inert fermions and scalars to realize the model, leading to two component dark 
matter candidates. 


The existence of new fermions and scalars provides rich phenomenon. Tiny neutrino masses are ex¬ 
plained with one-loop induced linear-seesaw-like mechanism. The charged-scalar mediates lepton flavor 
violation (LFV) of charged leptons. The relic density and the direct detection of the two component dark 
matter are investigated. The properties of discovered SM Higgs will be changed by the new particles. And 
these new particles provide plenty of new signatures at LHC. Especially, multi-lepton signals with missing 
transverse energy can be used to test our model. We find that our model can satisfy current constraints 
from phenomenons mentioned above. 

The rest of paper is organised as follows. In Sec. HU we introduce the realization of radiative linear 
seesaw and multi-component dark matter from gauged U(1)b-l- In Sec. Hill we discuss the phenomenon 
of lepton flavor violation, dark matter and collider signatures. Conclusions are given in Sec. HY1 
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TABLE I. New particles content: Gsm x U(1)b-l 


II. MODEL 

A. Model Setup 


In our model the neutrino masses are generated via the diagram depicted in Fig. Q] The new particles 
content and their charge assignment are listed in Table. Q] We add Asp generation Weyl fermions T Rr , T Rl , 
N\ right-handed Majorana neutrinos N Rn , and Afe pairs of right-handed Majorana neutrinos N' R and N” { to 
the SM where i. a and (3 are the generation indices. All the new fermions are singlets under SM gauge group. 
Five new scalars 771, S 2 , rj 2 , s 2 and a are also added to SM. Because the new fermions are all SM singlets, 
the B — L gauge symmetry satisfies all anomaly cancellations except for [U(1)b~l] x [Gravity] 2 and 
[U(1)b~l] 3 f! 1 Bn . Considering the conditions for the absence of [U(1)b~l\ x [Gravity] 2 and [U(1)b~l] 3 
anomaly, one has 


3 + (——)iVi + XN 2 + (—1 — x)N2 + (—1 )N q j ) — 0 (3) 

3 + {~) 3 Ni + x 3 N 2 + (-1 - x) 3 N 2 + (-1 ) 3 N^ = 0 (4) 

After solving the anomaly free condition, one obtains 

Ni = 2, N 2 = 1, = 1, x = ——— (5) 


Thus we have the some inert particles classified into two parts. In the first class, there are two Majorana 
right-handed neutrinos (Nr 1 , NR 2 ) and the inert scalars ( 771 , si) with their B-L charge being — In the 
second class, we obtain a pair of Majorana right-handed neutrinos (N' R , N R ) along with the inert scalars 
( 772 , .s 2 ) whose the B-L charges are irrational numbers or One notices that the new particles 

with both — i and the irrational numbers can not decay into SM particles. Therefore the lightest particles be¬ 
longing to the same class is stable and can be regarded as a dark matter candidate. The relevant Lagrangian 
for Yukawa sector is given by 

-Ly =yiLii’ R iT 2 §* + y'^L^RC + h a NR a i/j L si + f a iL^N Ra iT 2 rjl + ^Y Q N^ a N Ra a 
+ hN R il>BS 2 + fiL\N R iT 2 T \2 + —k N R N R a + h.c 


( 6 ) 

























5 


Without losing generality, we work in the basis where the mass term of Nr i j2 is diagonal. As for the mass 
term of N' R and N R , one can redefine the fields as 

Xi = 1 ={N'r + N^) X2 = A=( N ' R -Nr) ( 7 ) 

so that 

\YN%N' r <j \Y(&xx + rtX2)* (8) 

Now we get two Majorana neutrino eigenstates having the same masses. Note that there is no interplay 
Yukawa terms between Nri^ and (N R , N R ) because of the the B-L charges assignment they have. 

The scalar potential in our model is given by 

V(^,a,r)i,s 1 ,r]2,s 2 ) = - + A$(T> t T>) 2 - nl\a \ 2 + A ct |cj | 4 

+ v 2 m v \vi + \dvlm) 2 + tf 12 vlv2 + \ 2 (vlv 2 ) 2 
+ l s i| 2 + A Si | si | 1 + /r 2 2 |,s 2 | 2 + A S2 |s 2 | 4 + A SlS2 |si| 2 |s 2 | 2 
+ A^^-T) (77*771) + A^r/j^X^V) + A % $(T> t $)(r ? |? ?2 ) + A'^ (77*$) ($*772) 
+ 

+ As^lsil 2 ^^) + ^sirji l^i I 2 (^?i ^?i) + ^SlT]2 |si| 

+ A S2 d>|s 2 | 2 (^ t ^) + ^ S2T)l N 2 (7?lfil) + A S2 ^2 \s2\ 2 (vtm) 

+ A CT <j,|<T| 2 ( < f» t< f>) + \a m W\ 2 + A CTr?2 |a| 2 (?7^2) 

+ A si0 -|s 1 | 2 |cr| 2 + A S2cr |s 2 | 2 |cr| 2 + (fiisl&rjx + ^sl^V 2 + h.c ) 

(9) 

where // 2 . // 2 ;| , /./, 2 2 , //'^ and // 2 2 are taken as positive values and the value coupling constants fi\ and 
fi 2 in trilinear terms can be set as positive by re-phasing si and s 2 . Notice that there is no terms like sicr 2 
or s 2 (j 2 appearing in the scalar potential. This has two fold meanings: First, the inert scalars // 1 . 2 and si i2 
do not acquire the VEV after the SSB of <3? and < 7 ; Second, there exists a residual Z 2 x Z 2 symmetry under 
which all the inert particles are odd even after the breakdown of B-L symmetry. Therefore the residual 
Z 2 x Z' 2 symmetry stabilizes the inert particles, makes them to be two component dark matter candidates. 


B. Matrices of Scalar Particles 


After the SSB, the scalar $ and a is parameterized as 


T> = 


G+ 

v/2 


_ v a + cto + iGu 

■ n 


( 10 ) 
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where v$ — 246GeV is the VEV of the SM higgs doublet scalar and v a is responsible for the SSB of B-L 
symmetry [19]. The Nambu-Godstone bosons G + , G 0 and G a are absorbed by the longitudinal components 
of W, Z and Z' gauge bosons. For simplicity, we ignore the kinetic mixing between U(l)y and U( 1) b~l 
gauge boson Il20ll . Therefor the VEV v a provides a mass of U (1) b~l gauge boson Z' as Mz> = gs-LV a , 
where gs-L is the U(1)b~l gauge coupling constant. For the extra gauge boson Z', LEP-II provides a 
combined bound Mz> / gB-L > 7 TeV 12111 . which is just the lower bound on v n . Then we obtain the mass 
matrix for CP-even scalars (Ih) and op 


M 2 { 4 >q, op) = 


cos 9 sin 9 
— sin 6 cos 9 


where h stands for the SM-like Higgs 122 . 
masses respectively as 


Ml 0 

0 Mjj 


cos 0 — sm 0 
sin 6 cos 9 


(ID 


23] and H is an extra CP-even Higgs boson [241426] with the 


,2y2 


M'l — + X a vl — (Xv^ - X a v%) 2 + A^r^ 

MfJ = + ^CT V 1 + \J 

and mixing angle 9 determined as 


sin 29 = 


2Acr 

M 2 h ~ Ml 


( 12 ) 

(13) 

(14) 


On the other hand, the inert scalar (771, si) and (7/2, s 2 ) do not mix with <I> and cr due to the residual Z 2 
symmetry. The mass matrix for inert scalar fields are 


M{rt 1 ,si,rj2,s 2 ) = {v\, 4 , 4 ) 


^ Mu M \2 0 0 ^ 


M 

M21 M 22 0 0 


Si 

0 0 71^33 A/34 


m 

^ 0 0 A/43 A/44 j 


\S 2 J 


where 


+ 2^ lV 4 > + 2^ VlV ^ + 9 ^ 1 ^ 

M2 2 = + 2^ s i 4 > v \ + 2^ SiaV ° 

M33 = g 2 2 + + -X'^v^ + -Xl m v 2 a 

M 44 = + -Xs^v^ + -X S 2 C 7 i’l 

M \2 = M21 = 

V 2 

M34 = M43 = 


(15) 


( 16 ) 
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There is also no mixing between (771 , «i) and ( 772 , S 2 ), therefor a residual Z' 2 symmetry between the two 
classes can be realised. After diagonalizing the mass matrix, we obtain the mass eigenstates of inert scalars 


as 


A\ 2 \ / cos 0i, 2 — sin#i )2 

H h) \ sin cos^i^ 


rfl 
VI,2 


’ 1,2 


„• o a _ ^1,2^0 

Sin ^^12 r‘) f O 

’ Ml - M 2 

Jti 1,2 


(17) 


where 


m 1i, 2 = \ (K* + K, a + ^ M2 ^- M k,) 2 + 2 ^A) 


(18) 


M\ 


Hi ,2 



Here M m = Mu, M Sl = M 2 2 , M m = M 33 and M S2 = M 44 . 


(19) 


C. Neutrino Mass 
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FIG. 1. The one-loop diagrams for neutrino masses in our model. 


As shown in Fig. [H the tiny neutrino mass are generated by the linear seesaw mechanism except that 
the fiL terms are induced by a one-loop diagram. The effective mass matrix for active neutrinos depicted in 
Fig. Q]is expressed as 

Mull' = Ml u , + M^h, ( 20 ) 


where 

M lw 

M lu> 


v<p sin 0i cos 0i 
16vr 2 x/2 M^ 

v ^ sin 02 cos 02 
16tt 2 V2M 4 , 


Vl Y. hifu'Mj 


i= 1 

VihfvM x 


M A 2 


VMf - M\ x 


In (—-^r) — 
V M? ’ 


tj 2 Ml 1 

ln + 0 


m?-m 2 Hi 


Ma 


-.- A 2 Ml . M 2 AT 2 1 

spir 1" (4) - ^44 * (if)] +«« 0 


(21) 





















where M, (i = 1,2) denotes the masses for Nr\ and Nr 2 , M x denotes the masses for the eigenstates of 
N' r and N r . Tiny neutrino masses can be obtained using the following benchmark points: 


m = H 2 = 0.1 GeV,y = h = 0.0028, / = 0.01, = 300 GeV 

M Nri = 149.5 GeV,Mjv H2 = 200 GeV, M x = 150 GeV (22) 

M a o = 300 GeVjM^i = 270 GeV, M H o = 1000 GeV 
M a o = 700 GeV.M^t = 690 GeV, M H o = 62 GeV 

with the index of Yukawa couplings suppressed for simplicity. Then we get = 0.0164 eV(y / Am 2 3 ). 
The benchmark point given above seems rather unusual since the M v becomes a rank-1 matrix as 

(l 1 A 




D 


I 1 1 
1 1 1 


(23) 


which is obviously not consistent with the results of neutrino oscillation experiments. However, one reminds 
the expression of matrix (1231) is just the so-called flavor democratic model studied by many authors and 


related to some flavor symmetries [280. The matrix D can be diagonalized as 


V, t DV„ = 


( 


0 0 0 
0 0 0 
0 0 3 


\ 


f V3 1 V2\ 


= D 


Vu = 


V6 


-V3 1 y/2 

0 -2 y /2 ) 


(24) 


The unitary matrix V v corresponds to the democratic mixing pattern. Since I) contains a dominant non¬ 
zero element at (3,3) position, the flavor democratic structure in M u in eqn (l23l) can be viewed as a good 
approximation for our rank -2 neutrino mass matrix which exhibiting the strong and normal order of neutrino 
mass spectrum, i.e m± = 0 <C m 2 = \J A m'f 2 <C m 3 = \J A m'{ 3 . Thus the benchmark point we take is 
reasonable. It is noted in the flavor basis the democratic mixing matrix V v has already not been consistent 
with the PMNS matrix Upmns measured by experiments, however the V v can be corrected by the charged 


lepton sector V) to fit the neutrino oscillation data li2911 . 


The small values of pi 2 lead to sin 2#i 2 = (3.8, 7.2) x 10 -5 , which plays a key rule in the suppression 
of tiny neutrino masses. The choice of values of /i\ A is mainly for phenomenological consideration. First, 
in case of scalar DM, we have Z — S — S* coupling proportional to sin 2 The spin independent 
elastic cross section of DM requires sin #12 < 0.05130], setting an upper limit on /112 ~ (10 GeV) 
for electroweak(EW) scale inert scalars. Second, in our benchmark point with inert particles at the EW- 
scale and Yukawa couplings of the order 10' 2 or 10 -3 , rich phenomenon of new physics are expected for 
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LHC and LFV processes. On the other hand, larger values of // 1 ; 2 is possible if we decrease the Yukawa 
couplings to be more smaller values. But this will predict a too small branch ratios for LFV processes. Too 
small Yukawa coupling also seems unnatural from the viewpoint of model building. Another solution is 
to increase the mass of ^ or inert particle to TeV-scale, which is beyond the reach of LHC. For the DM 
candidate Nri^H®), its mass is set to be about half of the mass of the s-channel mediator If (h). Then one 
obtains the large enough DM annihilation cross section to account for the relic density. For the heavy dirac 
fermion if) and inert doublet scalars 771 , Aq, we choose their masses around 300 GeV, so that they are testable 
at LHC. The other inert scalars are around TeV-scale aiming to suppress the value of neutrino mass. 


III. PHENOMENOLOGY 


A. Lepton Flavor Violation 


The Yukawa interactions of charged-scalar rf 1 will contribute to the LFV processes of c 


Detail studies on LFV processes in scotogenic models [9] have been carried out in Ref. 113 ll. 


larged leptons. 


32]. Currently, 


the most severe constraint coming from MEG collaboration on muon radiative decay with an upper limit 
BR(/i —>• e-/)< 5.7 x 10 -13 (90% C.L.) [33]. In our model, the analytical branching ratio of // —>• ey is 


calculated as B ll. 


m 


BR(jU —>• ey) = 


3a e 


64itG 2 f 


E fififie 7-1 

~m\ f 

1=1 r]{ 



, Mi * p 
+ m\ F 

Vi 




(25) 


where the loop function F(x) is: 


F(x) = 


1 — 6x + 3x 2 + 2x 3 — 6 x 2 ln® 
6(1-a:) 4 ' 


(26) 


The benchmark point in Eqnl22l predicts BR(// —» ej) = 9.8 x 10 14 , which satisfies the cuiTent limit 


and is in the reach of future sensitivity 1331 . The limits on r observables are less stringent [ 361 


13711 . With 

the much natural Yukawa structure in our benchmark point, the predicted BR(r —>• irf) = 5.8 x 10 -13 is 
far beyond the future sensitivity. But on the other hand, hierarchal Yukawa structure f lt < \fi F \ < | f lT \ 
with f iT 0(1) is still allowed from phenomenological point of view. In this case, fermionic dark matter 


candidate F annihilation in the mass region between 2 GeV and 3 TeV through Lchanncl exchange of 7 / 
can satisfy dark matter relic density bound [32]- And as a consequence of hierarchal Yukawa structure, dark 


matter F annihilates mainly into third-family leptons: r + r and v T v T . 
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B. Dark Matter 


In our model, a multi-component dark matter scenario is possible due to the residual Z ' 2 symmetry 
between two sets of new scalars and fermions. For instance, we choose the lightest fermion Nr\ (refer as 
F) in set one and the lightest scalar H® (refer as S ) in set two as dark matter candidates. These dark matter 
candidates must satisfy two experimental constrains : (1) the dark matter relic density observed by Plank 


138] QomIi 2 = 0.1193 ± 0.0014; (2) cross section for direct detection of dark matter scattering off nucleon 


set by LUX 


H 


Theoretical calculation of dark matter relic density is well described in [39], and it is calculated with the 


help of packages Feynrules [40] and micrOMEGAs 1141 ll in our analysis. Because the i-channel Yukawa 


portal may suffer constrains from LFV or neutrino masses, and in our benchmark point, the f-channel con¬ 
tribution to relic density is less than 1%, thus we will focus on the s-channel h/H/Z' portal for simplicity. 
Firstly, the relic density of fermion/scalar dark matter for one dark matter candidate is presented in Fig. [2j 
where we neglect the conversion FFtt SS* between the two dark matter candidates. For the fermion dark 
matter F, the light h portal can not acquire a sufficient annihilation cross section, while the heavy H portal 
is still promising when Mp ~ Mh/% which is because the suppression of large v a = 8 TeV. Anyway, 
the Z' portal can easily satisfy the relic density when Mp ~ Mz> / 2. For the scalar dark matter S, it is 
dominantly made from the singlet scalar due to small mixing 62 . The relic density can easily be attained 
when Ms ~ M/,/2 and Ms ~ Mh /2, while the Z' portal is not promising, mainly because the small B — L 
charge of S and suppression of heavy Mz 1 ■ 

Secondly, we take into account the conversion of two component dark matter FF <->■ SS*, which can be 
mediated by s-channel h/H/Z ', where the //-portal is expected to be the dominant one. Therefor, HSS* 
and HFF are the two most relevant couplings to study conversion. For simplicity, we further assume 
A ar]2 = A S2(T = A, which determines HSS* coupling, and fix other parameters as discussed in FIG. [2]if not 


mentioned. 

The dependence of F/S relic density on A is shown in FIG. [3] For the fermion dark matter, when 
Mp > Ms, the larger A the larger FF SS* annihilation rate, and therefor the smaller the relic density. 
It is clear that the Qph 2 can differ by about one order of magnitude between A = —0.001 and A = —0.05. 
But when Mp < Ms, the effect of conversion FF —>• SS* is quite small. Noticing for A = —0.05, the 
H —> SS* can greatly enhance the total decay width of //, which causes the increase of Qph 2 around 
Mp ~ Mp/ 2. For the scalar dark matter, the increase of A will decrease the relic density significantly due 
to the increase of SS* annihilation. But, the friction of conversion SS* —>• FF keeps the same, since the 
Yukawa coupling HFF is fixed by Mp. The arguments are true if we exchange the roles of F and S, and 
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FIG. 2. Relic density of fermion (left) and scalar (right) dark matter as a function of Mdm for one dark matter 
candidate. Here, we set the relevant scalar interaction coupling X V 2 <s, = A/ 2$ = A S2 $ = —0.001, X ar)2 = X S2<7 = 
-0.001. We also fix M h = 125 GeV, M H = 300/600 GeV, sin 6 = 0.3, M z , = 4 TeV, and v a = 8 TeV. 


fix HSS* coupling while verify HFF coupling. 



MpfGeV) Af s (GeV) 

FIG. 3. The dependence of F/S relic density on A. 


Another aspect of conversion is the masses of the two dark matter. The left(right) of FIG. ED shows the 
F(S) relic density for Mg^ = 60,150,300 GeV with A = —0.02. In case of fermion dark matter, it is 
clear that the smaller the Mg the larger FF —» SS* annihilation rate, and therefore the smaller the relic 
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density. For relative heavy Ms = 150,300 GeV, the conversion has tiny effect on the F relic density when 
Mp < Mj{/2. In case of scalar dark matter, the dependence of 0.<,7; 2 on Mp is a little complicated, since 
the HFF coupling is directly related with the Mp. For Ms < 80 GeV, the effect of conversion is relatively 
small, and one expect that the larger the Mp the smaller the S relic density, which is mainly caused by the 
increase of H FF coupling. In medium mass region 80 < Ms < 200 GeV, the conversion effect would be 
dominant, thus the smaller the Mp the smaller the S relic density. In the high mass region, the conversion 
effect is comparable to HFF coupling effect, which makes the dependence of Og/i 2 on Mp nonlinear. In a 
word, the HSS* and HFF couplings play vital importance in dark matter conversion. The conversion can 
take place in both direction FF SS* and SS* —> FF when Mp ~ Ms, which can be obtained when 
both F and S are mainly annihilation through H -portal. If not the case, only the conversion of heavier one 
into light one is relevant B- 



M f ( GeV) 



M s ( GeV) 


FIG. 4. The effect of two component dark matter conversion on F/S relic density for fixed A = —0.02. 


Finally, we discuss the constrains from direct detection of dark matter. The Current experiment con¬ 
straints assume the existence of only one dark matter specie. However two-component dark matter candi¬ 
dates are predicted in our model. Therefore the contribution of cross section on nucleon for each specie 
should be rescaled by the fraction factor of relic density. We define the fraction of the mass density of z-th 


dark matter in case of multi-component dark matter (43, 


44fl: 


0,7; 2 


G; = 


O CDMh 2 ’ 


(27) 
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where i = F,S in our consideration. Thereafter, the up limit of direct detection is: 


£F 


es 


cr, 


exp 


i, T -&F—N H T~j—&S—N < ,, 

Mf Mg Mom 


(28) 


Here, <7F_jv(<TS'-iv) denotes the scattering cross section of F(S') with a nucleon N. The benchmark point 
in Eqn. [22] gives the spin-independent scattering cross section = 1.10 x 10 -46 cm 2 = 

1.62 x 10” 44 cm 2 ) with VLph 2 = 1.12 x 10~ 4 (flg/i 2 = 2.64 x 10~ 4 ). Although the bare Og-N l ar g er 
than the LUX upper constraint 1.1 x 10 _47 cm 2 /GeV[3(j], the contribution of scalar S to the scattering on 
nucleon is suppressed because of its small faction eg = 2.21 x 10~ 3 . The value of expression on left of Eqn. 
[281 is 1.3 x 10 48 cm 2 /GeV, which is smaller than current LUX bound. Thus, the fermion dark matter is 
dominant in this scenario, while the scalar dark matter must be less than 4% to escape current LUX bound. 

C. Collider Signatures 


As shown in our bench mark point (Eqn. 1221) . the new particles are all at the electroweak scale, which 
makes them testable at LHC. Interactions between these new particles and SM Higgs h will of course 
modify the properties of h, thus give us some indirect hints. Nowadays the most precise measurement of 


M/ ( is the combined results of ATLAS and CMS [45]]: 


M h = 125.09 ± 0.21 (stat.) ± 0.11(sys.) GeV 


(29) 


Apparently, the extra of new scalars and fermions would change the decay rates of SM Higgs h. Lor in¬ 
stance, mixing between h and additional scalar singlet H will modify tree-level h decays. And the additional 


charged-scalars rjf 2 will contribute to the loop-induced decays as h 77 [46]. 


It’s well known that, for Higgs-portal dark matter, upper limits on Higgs invisible decay is also inter¬ 


preted an upper limits on dark matter-nucleon scattering cross section 14711 . Direct measurement of Higgs 


invisible decay in associated with Z by ATLAS set an upper limit of 75% at 95% C.LM48]. Combined 


analysis with Higgs signal strength give a more tight upper limit of 37% at 95% C.L. [49 


51]. In the future, 


the weak boson fusion channel might have the ability to probe invisible decay to 2-3% with 3000 fb _1 at 
LHC 1152ll . Light dark matter candidates in our model will contribute the the invisible decays of h in our 
model. In case of scalar dark matter, the branching ratio of Higgs invisible decay, i.e., BR(/i —>• SS*) is 
1.7% for Ms = 62 GeV with A S2 $ = A S2(7 = —0.001, sin 9 = 0.3. On the other hand for fermion dark 
matter, BR (h —> FF) is 1.7 x 10 -6 for Mp = 62 GeV with v„ = 8 TeV, sin 9 = 0.3, although it might 
not be favored by the constrains on relic density of dark matter. 


Then we discuss tl 


0.23 at 95% C.L. |53 


ic mixing between h and H. The analysis of signal strength of h constrains sin 2 9 < 


54 ]. Direct search of H in the ZZ and WW channel now has push this limit down to 
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FIG. 5. Branching ratios of H as a function of Mr for sin 0 = 0.3(left) and sin 6 = O.Ol(right) in our benchmark 
point in Eqn. [22] 


sin 2 9 < 0.1 with no new physics contribute to the decays of H [55]. Future hadron collider, i.e., HL-LHC, 
has the ability to probe sin 2 9 ~ 4 x 10 -2 , and lepton collider, i.e., CEPC, could reach sin 2 0 ~ 2 x 10~ 3 
[56]. In Fig. [5] we show the branching ratios of H for two values of sin 0(0.3, 0.01). For a relatively 
large mixing angle sin 9 = 0.3, the heavy neutral Fliggs H decays dominantly into SM particles. The 
branching ratio of invisible decay H —>■ SS* can reach 10% for Mu ~ 165 GeV, and we expect it becomes 
dominant when Mu < 160 GeV. The branching ratios of H decaying into other new physical particles are 
below 10~ 3 in this case. While for Mu Myp, it is well known that decays of H into vector bosons are 
determined by their Goldstone nature, which implies: 


BR (H -»• hh) » BR (H -> ZZ) « ^BR (H -»• WW) 


(30) 


The asymptotic behavior of this relation is clear shown on left picture of Fig. [5] On the other aspect, for a 
relatively small mixing angle sin 0 = 0.01, decays of H into SM particles will be suppressed and decays 
into new particles will be greatly enhanced. H —>• SS* is dominant when Mu < 400 GeV. The branching 
ratio of will reach about 0.25 when Mu ~ 600 GeV, which is comparable with H —> W + W~. In this 
case, H FF is below 10% and H —> r] + rf~ / A^A®* is below 2%. 

The heavy neutral Higgs H is testable for large mixing angle 9. For example, the promising channels to 
probe heavy neutral Higgs H would be ZZ — > 4 1, ZZ — > 2l2v, ZZ — > 2l2j, ZZ —> 212t, WW M 2l2u, 
WW —> lv2j, hh M 4b, and hh M 2h2^ js?]. At the same time, we would like to mention that the heavy 
Higgs H could enhance the di-Higgs production hh [58] by a factor of 18 comparing to standard model 
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case [590. For small mixing angle 9, production of H will be suppressed by this small 9, thus makes it 


challenging to probe directly at colliders. 


CU 

O 



FIG. 6. Total decay width of Z' as a function of M z' (for fixed values of fjri-L), and gs-L (for fixed values of Mz*)- 


Next we review the properties of U(1)b~l gauge boson Z' . With about 20 fb 1 data at 8 TeV LHC, 
bound on Z' has been push up to 2.95 TeV by CMS through the ratio R a = a(pp —» Z' —> £ + £~)/a(pp —> 
Z —>■ £ + £~), where l = e, p B60h - In our benchmark point, we choose Mz> = 4 TeV and gs-L = 0.5 



In table [Til we give the decay branching ratios of Z' in our benchmark point. The dominant decay 
channels of Z' are qq, ll, and while all of the new particle final states only account for about 20%. A 
distinct feature of (7(1) b~l gauge boson Z' is the definite relation between quark and lepton final states: 


BR(Z' -)• qq) : BR(Z' ll) ~ 2 : 3 


(31) 


after summing over all flavors. This relation can be used to distinguish f7(l)^-L gauge boson Z' from 
Z' in other models fl- More practical on experiment, the B-L nature of Z' can be tested if BR(Z' —> 
bb)/BR(Z' —>• [i + fi~) = 1/3 is confirmed [13]. In our model with only left-handed light neutrinos, the 
dominant invisible decay channel of Z' is BR(Z' —>• VfVr,), which is half of BR(Z' —>• ll). Further with 
dark matter candidate in our model, Z' invisible decays get additional contributions from Z' into dark 
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matter pairs. For instance, BR (Z 1 —> inv.) could be 0.2457, 0.1990, 0.1964 for FF, FS and SS dark 
matter separately. So precise measurement of BR(Z' —> inv.) will shed light on the nature of dark matter. 


qq 

ll 

VLVL 

ipip 

NriNri 

Nr 2 NR 2 

Xi Xi 

X2X2 

0.25 

0.38 

0.19 

0.063 

0.0077 

0.0077 

0.016 

0.032 

HH 

hh 

AlA"* 

H\H\* 

too 

too 

* 

H°H°* 

ViVi 

vt % 

0.030 

0 

0.0076 

0.0051 

0.0010 

0.0013 

0.0076 

0.0010 


TABLE II. Branching ratios of Z' in our benchmark point. Here, we set Mu = 300 Ge V and sin 9 = 0 for simplicity. 



FIG. 7. Branching ratios of ip as a function of M^. We have fixed sin0 = 0.3, Mh = 300 GeV and y = h = 0.01. 


Another interesting feature of our model is the existence of heavy Dirac fermion ip, thus there are no 
lepton number violation (LNV) decays as ip ^ W~l + . For AT,/, < My, the Higgs decay into a pair of 
light and heavy neutrinos, h vip + ipv will open, which could increase F/, by up to almost 30% and 


significantly affects Higgs searches at the LHC [6fD- In this paper, we consider AT,,, > Mh. Therefor, 
decay channels of ip could be W + l~, Zv, hv, and if kinematically allowed Hu, A®Nj, A®Xj, TT V °Aj, 
= 1) 2) are also possible. Due to tiny mixing angle 0\ 2 , branching ratios of ip —>• A®Nj, A®Xj 
are negligible. In Fig. [7J we show the branching ratios of ip. It is clear that ip will decays dominantly into 
standard model final states for comparable Yukawa couplings of y and h. Approximately for AT,/, A> M\y , 
we have: 

—^tj—BRI^ —»■ hu) « BR (ip —>■ Zv) ~ -BRf^ —>■ W + l~) 
cos^ 0 2 


( 32 ) 
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Decays of ip into new physical particles is small in this case. BR(?/> —>• H^x) is about 10%, once is 
kinematically opened. BR(-0 —>• Hu) is suppressed by sin 2 6 , thus it is always much smaller. As shown 
in table |TTJ BR (Z' —>• ipip) — 0.063 for one generation in our model, so ipip can be produced through Z’ 


portal. A possible promising signature is the tri-lepton channel j62ll : 

pp ^ Z / ^ W + l~ + W~l + -> 2 iH^jj + Xt 


(33) 


The cross section of this tri-lepton signal is about 0.017 fb in our benchmark point, so the tri-lepton is only 
promising at future high-luminosity LHC. The mass of ip can be reconstructed using the transverse mass 


of two opposite sign leptons with missing transverse momentum [62]. Another feature of ip is the possible 
large mixing with ui comparing to canonical type-I seesaw |3j]. As discussed in Sec. Ill C[ the mixing V vy) 
between vl and ipr is Md/^/Mq + M^. For Mr, ~ 0(l)GeV, ~ 0(lOO)GeV, V m ^ ~ 0(1CT 2 ). 
Thereafter, ip could be largely associated produced with charged leptons through W [64]: 


pp —>■ W* —»• l±ip -> / ± + W ± H -> l^Fjj 
pp — W* —»• l±ip -> + W ± P -> 


(34) 

(35) 


The production cross section cr(l±ip) = 350 x Vj y . \ 1 fb in our benchmark point. And it might be promising 
at 14TeV LHC with about 100 lb 1 . Testability of this heavy Dirac neutrino ip is less promising than the 


68Q. 


heavy Majorana neutrinos with same mixing scale, since the latter could rise LNV signatures 165 - 

Finally, we discuss the decays of inert scalars and fermions. Nr i and H® are a dark matter candidate 
as in our benchmark point in Eqn. [22 ] Decays of Nr 2 are dominated by Nr 2 —> —>• / ± / 1 Nr\ and 

Nr2 — >• uA'l* — y vvNr i through the Yukawa coupling /. Decays of Xi are Xi H^ip* with the off-shell 
ip* further decaying into W + l~ /Zu/hu/H^u. A 1 ) and pf mainly decay through the Yukawa coupling /, 
which leads to A® —>• uN[,, and pf —>• I^Nr^. The heavy 2L odd scalar H® decays into ipNjn through 
Yukawa coupling h a and into h A ( ) through trilinear coupling p\ . Similar, decays of A ( ) are ->• lJ Xi and 
A° —>• hH ®, while decays of pf are pf —>• l±Xi and vt 


Particles 

vtv i 

vt A ° 

A°A°* 

vtv 2 

vt A °2 

a0 a0* 
^2^2 

a (in fb) 

5.8 

16 

3.6 

0.089 

0.31 

0.075 


TABLE III. Production cross sections for inert scalar doublets. 


The inert scalar doublets can be pair production through DY process. In Table [Till we list the production 
cross sections for inert scalar doublets. Many signatures can be risen from the inert particles. In Ma’s 
scotogenic model |^[, promising signals of the doublet scalar on collider are multi-lepton final states with 
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missing transverse energy 1169 


7ill . Similar signals can also be produced in our model, for example: 


21 + Xt '■ ViVi l + ^m + l Nri (36) 

: rfirii -> W + H%* + W~H$ -> l + viH%* + l~viH% 

31 + Xr : vtA -> / + iVi? 2 + i/JVjn -> I+I^FNri + i/JV m (37) 

: jjf -> + h-H® -> / ± ^,77 2 ° ( * ) + 

4 1 +Xt ■ vtvl ->• / + JV/a + riVflt ->• (+/ ± ( =F iV R1 + 1~Nri (38) 

: -)• hH^H^ -> l + ri + l~H% + H$* 


With much different decay topologies between our and Ma’s model, it would be distinguishable even with 
same signals. Apart from these multi-lepton signals, there are also some other interesting signals in our 
model, i.e.: 

2 l±jj +Xr : rgA% ->• + /ttf 2 0 -> l^H 0 ^ + l±v t jjH\ 2 ° (39) 

/ ± 66 +Xt : -> W ± H^* ) + hH$ -»■ / ± ^77 2 ° ( * ) + bbH$ (40) 

The lepton number violation signal 2/ ± jj + Xt suffers much lower SM background, thus might make this 
signal very promising on LHC. The ^ftfr+Xr has a relatively large production rate due to h —>• bb dominant 
in h decay, so it might also be promising. 
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IV. CONCLUSIONS 

In usual canonical seesaw mechanisms, it requires the heavy states with the scale of masses being grand 
unification scale to generated the small neutrino mass. In linear seesaw scenario with m u — n lMjj/M q,, 
the neutrino masses suffers a two fold suppression by both lepton number symmetry violating term ///_ and 
heavy mass Afq,. The linear seesaw model can lower the seesaw scale such that new physics may arise at 
TeV scale. In this work, we construct a radiated linear seesaw model where the naturally small term /i£ are 
generated at one-loop level and its soft-breaking of lepton number symmetry attributes to the SSB of B-L 
symmetry at TeV scale. To satisfy the anomalies cancelation, the value of B-L charges for inert particles 
are found to be exotic such that there exists residual Z 2 x Z' 2 symmetry even after SSB of B — L gauge 
symmetry. It is shown that the residual symmetry stabilizes the inert particles as dark matter candidates. In 
our model, we introduce two no-interplay classes of inert particles to realize the model such that the lightest 
inert particles belonging in each class play is the dark matter matter candidate. Therefore we have propose 
a two-component dark matter model. The seesaw scale of radiated linear seesaw scale can be as low as a 
few hundred GeV, leading to interesting phenomenology. 

Given a benchmark point at electro-weak scale, we illustrate the main prediction of our model. For 
the Yukawa coupling fu(fi) at 0.01 order, our benchmark point predicts Br(/i —>• e^y) ~ 10“ 13 , an order 
slightly under the current constraints and in the reach of the forthcoming experiment. The two component 
dark matter candidates are realized in our model. To account for the observed relic density, the annihilation 
of dark matter are dominant by the s-channel scalars h/H or gauge boson Z'. For the fermion DM, we find 
that the /i-channel is excluded. But it is still allowed for the //-channel and Z'- channel. On the contrary, for 
the scalar DM, the Z'-channel is excluded while the h/H -channel is allowed. And the heavy Higgs H also 
plays a vital import rule in the conversion between fermion and scalar dark matter. Collider signatures of our 
model are also very rich. The precise measurements of SM Higgs h will put tight constrain light scalar DM 
and heavy scalar H. With a relatively large mixing angle sin 0 = 0.3, the H —>• ZZ. W + W~ channels are 
testable at LHC. For the extra Z boson and heavy lepton ip, the tri-lepton channel of Z —>■ ip ip is promising 
at HL-LHC. With a larger cross section, the associated production of / ± L may be more promising. The 
inert doublet scalar can also produce multi-lepton channels. And some distinct channels, as 2£ ± jj + Xt> 
iHb + Xr, can be used to distinguish our model. 

Finally, we would like to mention that the radiated linear seesaw model we proposed is the minimal 
version where only one 'P fermion mediator is included. In this scenario, M v is a rank-2 mass matrix 
and the lightest neutrino must be massless. However, more complicated scenarios exist, corresponding 
to other solutions of the anomaly free condition. Then one may obtain the rank-3 neutrino mass matrix. 
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Such scenarios predicts more new particles with different B — L charges, the model construction and the 
phenomenology deserve us further study. 
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